SDF - INTRODUCTION 


A single degree of freedom system consisting of a 
mass mand a spring with stiffness k is considered. 


L» u(t) 
k k (N/m) 


pm p(t) 


The spring is undeformed for u = 0 


The equation of motion can be derived in two ways. 


Newton’s second law 


u(t) >0 u(t) 
—> 
k p(t) ku(t) p(t) 
pie TP 
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D’Alembert’s principle 


The system is supposed in dynamic equilibrium. 
The principles of statics are applied by 
introducing a fictitious inertia force, a force equal 
to the product of mass time its acceleration and 
acting in a direction opposite to the acceleration. 


Free-body diagram 


u(t) 
ku(t) =a 


oe 


mut) 
“fictitious inertia force 


Dynamic equilibrium 


p(t) — ku(t) — mi(t)=0 


Stiffness k 


The stiffness k has the same definition as in 
the displacement method: 


The stiffness k is the external force that is 
needed to keep the system in equilibrium 
when a unit displacement u = 1 is applied. 


et 


k 
abs 


spring external force 


Combination of two springs — case 1 


sed 


The stiffness for this 
system Is trivial: 
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Combination of two springs — case 2 


ic ae 


a 
Or ©) 


This case is more 
complicated. 


aa 
eae 
statics: R=hj=k 
d;+d> =1 aa eee 
k, ko 


Structural dynamics 1.3 
Some structures can be idealised as SDF 





Example aie S <6 The mass of this SDF system is m, the 
oi —> mass of the roof. 
4 
In statics, this frame 
has 6 active degrees 
Grirecdonn The stiffness is determined in the classical way: 
u=1 
2 : me 
Aa o k 
1 
By neglecting the 
axial deformations, 
3 d.o.f. disappear. 
1 
rigid beam 12et nA ; SEI 
| ig | 2 
—» | 
Only one d.o.f. is left if : 
the frame is consisting of 
an heavy roof supported pe saake = ( ae — 
by light columns. L 


Rotations 





massless 





A SDF system can ~ 
also have a rotational 
movement. 


1Td*G 


iInsialics: = M—Ce. c — 





eZ 


indynamics: | J6 + C@ = M(t) 


J is the moment of inertia. 


2 
foracircular disk: J= a 
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J@ is the moment created 


by the fictitious inertia forces. 





Example: Calculation of J for a bar 


. _——._ inertiaforce f, 
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SDF - FREE VIBRATION 


UNDAMPED FREE VIBRATION 
The structure is disturbed 
from its static equilibrium k 


and then vibrates without 
any applied forces. 


The equation of motion is: 


The solutionis: u(t) = Acos(@,t) + Bsin(a@,t) 


k/m| (rad/s) natural circular frequency 


A and B are determined by the initial conditions 


Ur_g = Up — Uo =A 


Structural dynamics 1.5 


u(t) = u, cos(@,f) + “0 sin(@,f) 
Op 


which can be written as_ | u(t) =Csin(@,t + 6) 


C— Ger (ey can cose = Lol n singe 


DB AEA. 
LY 


natural frequency 


© 


natural period 


27 l1 @ 
ar (s) ae OR (Hz) 
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Energy in undamped free vibration 


At each instant of time, the total energy E 
is made of two parts, the kinetic energy E, 
and the strain energy Es. 





E(t) = E,(t) + Es(t) 


Remark : the conservation of energy can be 
used to derive the differential equation. 


_ =m [C a, cos(@,t + 6)) + sk [Csin(@,t+ 6)/ 


= me? @; COS” (@pt + 4) ier C’ sin’(a,t + 0) ' ' 
2 2 E(t)=—mu?+—ku? 
1 2 2 
Ge (k= ma; ) dE 
2 conservation of energy —-> Ta 0 


oS ee T=0 mu+ku=0 
conservation of energy. AP uu=0 — mu+ku= 


F(t) is constant, which implies 





VISCOUSLY DAMPED FREE VIBRATION 


Friction in the structure is u(t) 

idealized by a linear viscous md 

damper which develops a c 

force proportional to the k me 

velocity =o 
f, =—cu(t) 


The equation of motion is: 


if C<c, =2Vkm (criticaldamping) the solutionis 


u(t) =e >" [Acos(a,t) +Bsin(a,0)| 


damping ratio 


damped pulsation 


Op = On 1-€ 
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A and B are determined by the initial conditions. 


u(t) = -Ev,e °"'[Acos(a,t) +Bsin(a,)| 


+e oun! [— Aa, sin(a,t) +Ba, cos(a,t)| 


U9 =U 2 U=A 


u(t e Sent 4 cos(a,t) + Ho 52h sino.) 


QD 





The solution can also be written as 


Uo 


Oy Co sin(@,t + 8) sind = a 





, Z 3 
C= +o sent) cose = Ho FSP at 


Wp Wp 





Decay of motion 


A free vibration test can be used to determine 
experimentally the natural frequency and the 
damping of a structure. 
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p periods between two maximal 
points u, and Un,» 


th+p =t,+ pl, 


ie Ce S°n'n sin(ao,t, +) 
Unrp Ce 2%nlntPld) sinfo, (t,+ pl>)+9) 
= ec@nP Tp 
27 


U 
In" =¢@,pT, =¢@,P 3 
Un+p O,\1-€ 





COULOMB-DAMPED FREE VIBRATION 


Coulomb damping = u(t) 


results from friction k Aaa 
; fae ry friction 
against sliding of HWS 


two dry surfaces. 


The friction force is F = uw N where uw denotes 
the coefficients of static and kinetic friction, 
taken to be equal, and N the normal force 
between the sliding surfaces. 


F is assumed to be independent of the velocity 
of the motion and its direction opposes motion. 


The equations of motion from left to right are 


ee 





u(t) = A, cos(@,t) + B, sin(@,,t) — F/k 


The equations of motion from right to left are 


mu+ku=F 


u(t) = A, cos(q,t) +B, sin(@,t) + F/k 
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The constants A,,B,,A,,B, depend on the initial 
conditions of each successive half-cycle motion. 


The plot of the solution is 





DAMPING 


Different damping models can be used, a 
viscous proportional damping is the most 
used approach. 


There are two reasons for that: 


¢ The mathematical equation which describes 
the motion is easy. 


¢ This model gives results which are often in 
very good agreement with experiments. 


A consequence is that the damping 


coefficient € can only be determined 
by experiments. 
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MULTIPLE DEGREES OF FREEDOM 


SDF mu+cu+ku= p(t) 


MDF = [m]{u}+[e]{u}+[k]{us={p(d)} 


[k] stiffnesss matrix 
[c] damping matrix 


[m] mass matrix 


The mass and _ stiffness matrices are 
obtained by (finite element) discretisation of 
the structure. 


The damping matrix cannot be obtained by 
discretisation, a different approach must be 
used. 





SDF - HARMONIC EXCITATIONS 


A harmonic load b> u(t) 
is applied to the 


k 
m 
structure. i 


Without damping, the equation of motion is 


mu+ku=p,sin(at) 


Example 1 


The system consisting 
of the mass m and the 
excentric mass M is 
considered for writing 
Newton’s equation. 


k/2 k/2 


(m-m)uU + m(V+U) = -ku 


mi + ku = Mra sin(at) 


Po sin(at) 
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Example 2 


y=U+Ug 


my=-ku 





= Ug(t) = Ugo sin(at) 


If uis to be studied (e.g. earthquake) 
m(u + Ug) =-kKu 


“ On 
> | mU+ku = Mugs @ sin(at) 


If yis to be studied (e.g. floor isolation) 


my =-—K( y—Ug) 


>} mytky = kuggsin(at) 
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WITHOUT DAMPING 
u(t) C and @are determined by the initial conditions 


L-» 
k Pp sin(at) : 
u(t) is a summation of two sinus and is not 
ee ae . mM 
MU+ K U= Dy SiN(@t) defined for @ = @. 


In reality, the damping implies that u,(?) 
disappears after some time. Then the solution 
( steady state response ) is 


The solution u(t) of the differential equation is 
the sum of two parts u, (f) and u, (0). 


u(t) = Up(t) + Up(t) pik 
homogeneous solution wu,(t) = Csin(a,t + 4) 


conclusions 
particular solution u,(t) = Asin(at) 


«After some time, the structure vibrates with 
the same frequency as the applied force. 


i,(t) =-Aw” sin(at) 


ay ~mA@?+kA= Pp, Ae Po = Po/k «The amplitudes of the vibration are infinite 
k-mo~ 1-(a/o,) when @ = @,. 





u,(t) is the transient response 


ee Po/k 
u(t) = Csin(@,t + 8) + 1—(/@,) u,(t) is the steady state response 





WITH VISCOUS DAMPING 


mu+cut+ku=p,sin(at) 


The homogeneous solution u,(f) (transient 
response) disappears after some time. 


u,(t)=Ce >" sin(@,t + 8) 


The particular solution u,(t) (steady state 
response) is of the form: 


u(t) = Asin(at — ¢) 
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after calculations (see the book), it is obtained 





p(t) = Poi k sin(ot—¢) 
-(o/0,? |’ +bélo/a,)]” 


tan = ee 0<6<180° 





Remarks 


* The total response is u(t) = u(t) + u(t). But 
after some time u,(t) disappears and u(t) = u(t) 
(steady state response). 


¢ After some time, the structure vibrates with the 
same frequency as the applied force. 





This numerical example shows that the 
transient response u,(f) disappears after 
some time and that only the steady state 
response u,(1) is then left. 


u(t) = Up{t) + Up(t) 


u(t) =Ce°™' sin(a,t + 6) 
Po Kk 
ll -(o/a,)]’ + Ré(o/o,)]? 


+ 





sin(at — @) 


C and @ are determined by the initial conditions. 


Structural dynamics 2.4 





total response u(t) 


Sen ae transient response u,(t) 


For this case, the steady state is obtained 
after about 20 seconds. 


Dynamic factor 


After some while, the structure vibrates with the 
same frequency as the applied force. It is the 
steady state response u,(1). 


The amplitude of these vibrations are now studied. 


Dy /K 





u,(t) = sin(at — ¢) 
p 7]2 2 
hi -(a/o,)7|° + 2E(o/o,)] 
The static deformation 220s 
due to a static load p, is Uso = he 


The amplitude of the vibration is equal to the 
product of the static deformation times a 





dimensionless dynamic factor Ry. 


Fq(5,@/@p) 
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R, can be plotted as function of the ratio w/a, 
for different values of the damping coefficient €. 


4 ea é=01 
3.5 
3 

Hae € =02 
2 

1.5 ¢=04 





(@/@,) 


o/0,<0.25 — R,=1 "quasi static"response 


Ari +eti+ ku=p,sin(at) > u=p,/ksin(at) 


@—o, the amplitudes of vibrations 
become large :Resonance 


Resonance 
Resonance is reached for 
ee 2 
O = 0, = Op, 1-26 


For this value of w, the dynamic factor is 


1 
if ¢<0.1 then Pamax ~ 52 and 


o, =/k/M = @,=@,\1-€" = @, = @,V1—2€" 





Example 
Gayo — leis = 25 


The deformations are 25 times the static ones. 
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Band-width method 
( Experimental method to determine ¢) 


The structure is excited by a harmonic load. The 
frequency of the load is increased step by step. 
At each step, the amplitudes of vibrations of the 
steady state response are measured. This 
implies that at each step, some time must be 
waited so that the transient response 
disappears. The curve A, as function of wis then 
obtained experimentally. 





Force transmission and vibration isolation 


A harmonic load is applied to a structure. This 
structure is connected to the ground through a 
support modelled by a spring k and a damper c. 


p(t) = py sin(at) 






U 
J Steady state response 


u(t) = = Rg sin(at — ¢) 


The force transmitted to the ground is 


f(t) = ku(t) + cu(t) 


= poRy sin(at — 6) + POE Rywcos(ot -¢) 
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f(t) = PoRalsin(at — ¢) + 2¢(@/@,) cos(at — 9)| 


= PoRq1+[2é(@/a,)]° sin(at- 9 + @) 


The transmissibility TR is defined as the ratio 
between the amplitude of the transmitted force f- 
and the amplitude of force applied to the structure. 


es 2) 


Po \[i-(w/a,)?|" + 2e(@/@,)]° 





TR is dimensionless. 


The objective is to choose the support (k, c) 


such that 7A is as small as possible. 





The mass m of the structure and the frequency 
of the load @ are imposed. The problem is to 
choose a support with k and &€ such that the 
transmissibility TR is minimal. 


TR as function of the ratio w/q@, is plotted for 
different values of ¢. 
Ok in 


The transmitted force is less than the applied 
one if the natural frequency @, Is such that 
O!10,>1.4. 


A low TRis obtained for low values of @, and ¢. 


However, a very low @, implies a low k and 
therefore a too large static displacement p,/ k. 


Besides, a very low €é& implies high 
displacement amplitude while passing through 
resonance which may occur before the load 
reaches the circular frequency @. 


Natural rubber is a good compromise and is 
often used for the isolation of vibrations. 
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SDF - ARBITRARY EXCITATIONS 


This chapter studies the response of a SDF 
system to pulse, impulse and periodic excitations. 


RECTANGULAR PULSE FORCE 


p(t) 


gop me p 
m p(t) 


initial conditions: U, =U, =0 


An undamped SDF system is loaded by a 
rectangular pulse force. For this case, the 
differential equation can be solved analytically. 
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NSPE, mu+ku=D), 
u(t) = u,(t) + u(t) 
u,(t) = Acos(qa,t) + Bsin(@,f) 


u(t) = Acos(@,t) + Bsin(@,t) + a 


Uy =0 A+ Po =0 
= k 


toe) Ba, =0 


u(t) = = [1 -cos(a,t)| 


{2 mu+ku=0 — freevibration 


u(t) = u(t,)cos(apt,) + wo) sin( cot) 


n 


en al ~co0s(@pt»)| 


U(ta) = PEEP sin( nto) 


u(t) = a [i - cos(@,t,)|cos(@,t — Opto) 


4 Fo sina) SIN(@pt — Opto) 
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u(t) = cs [cos(a,t — @pt) — COS(@pt,) COS(@,t — Opto) 


+ sin(@,t,) sin(@,pt — @pto)| 


by using cosacosb-sinasinb=cos(a+ b) 


it is obtained 


u(t) = = [cos(a,t — @,t,) — cos(a,t)| 


IMPULSE LOAD 


eee 
+» u(t) 
k 
pC) 
() () t 
iG 


An undamped SDF system is loaded by a short 
impulse force. An approximative response is to 
be calculated. 


du(t) 
dt 





Newton’s law for f<t, =m +k u(t) = p(t) 


Hypothesis: ft, is so small that the displacement 
uis still zero at t,. (U,= 0). 


t, is very small 
| . ! > u(t,)=0 


uo =0 


Then ku(f) can be neglected in Newton’s equation 
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lo 
du(t) = — p(t)dt eae == fp(i)dt 
0 


ii lo 
I = [p(t)dt 
0 


At t,, the impulse load is assumed to have 
produced an initial velocity, but no displacement. 


Newton's law for t2t, (free vibration) 


u(t) = u(t,) cos[o,(t—t,)] + wee) sinloon(t- to) 


. iS 7 I 
— the maximum displacement is | Ura, =—— 
MO, 


This result is only valid if tf is small enough so 
that u(t.) = 0. In practice it means ft, < T,,/10 


In such case, the maximal deformation does not 
depend on the form of the impulse load, but only 
on the value of the integral 7 
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PERIODIC EXCITATION 
Newton’s equation is 


ASDF system is excited fet “ a , ae 
aU | ae 2 EF COS(/@ot) + b; sin(ja,t)| 


by a periodic (but not m p(t) 
jel 
harmonic) load. ees 
The steady state response is calculated by 
using the theorem of superposition. 
Ao 


mu+ku=a — Up= k 





bi 
eee oo) 


iL 2 37, 
@ e & mu + ku = b;Sin( j@pt) > Up = ; 
ae 1—(j@_/ @p,) 
Idea : a periodic function can be separated into b./k 
mi+ ku = b;sin( j@,t) > u, =———’———- sin( ja, t 
J (J O ) p 1- (ja, /@,) (J Oo ) 


its harmonic components using Fourier series. 


aes la; COS( jot) + b; sin( ja,t)| 
fel 
The total steady state response is then 


p(t)cos(ja,t) dt 


1/k ee 
5 la; COS(/@,t) + 6; sin(ja,‘)| 


Dera) 
j=1 1—(J@g/@p,) 





To Dae ae 
J plt}dt — by= = J plt)sin( jargt) dt 
0 @ @ 


Example = o<t<7,  p(t)=2t 
i 

Ie anemia 

a=— | tdt=-2 

ane 

T, 
aj=— f Po tcos(ja,t)dt =0 
I, 0 Tp 


2 1a p ere p 
ae | ee are 


The total steady state response is then 


Da feds ae 
+ _sin( ja, ft) 
2k j= 1-(f@/@,)° ° 


u(t) = Po — 





“+p Do _ Po _P 
Wilh == 2 if) =O Se 
tt) 2k «1k 





2 ot) 





H sin(3@,f) i sin(4a@,f) 
—105 =)597 


 Snice! S sin(2@,f) 
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Only the first three Fourier terms in u(t) must be 
considered to get an error less than 2 %. The 
plot of the load approximations with 2 and 3 
Fourier terms shows that an_ inaccurate 
approximation of the load gives an accurate 
approximation of the response. The reason is 
that the higher frequencies in the load do not 
give any contributions to the response. 





NON PERIODIC EXCITATION 


p(t) u(t) 


t t 


Fourier series can also be used for non 
periodic loads. As example, the response u(t) 
of a damped SDF system loaded by the force 
p(t) is given above. 


p(t) and therefore u(t) are non periodic. They 
become periodic by artificially adding the same 
model. Then a Fourier analysis can be done. 


The Fourier analysis of the artificial periodic 
problem is performed by only considering the 
steady state response. But the result obtained 
for the real problem includes both the transient 
and steady state response. 
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Remark: u(t) must be zero at the end of the 
period 7,. 


Such an approach is not used for real problems, 
but it constitutes the basis for signal analysis, 


see the chapter experimental dynamics. 





RAYLEIGH’S METHOD 


Rayleigh’s method can be used to calculate 
approximately the lowest natural frequency of 
beams. 


INTRODUCTION 


u(t) 
Free vibrations of an k i 


undamped SDF system 
(o) en () 


are considered. 
u(t) = u, sin(@,,t) u(t) = u,@, sin(@,f) 


strain energy 


E.(t) = = Kutt)? 7 = kue sin? (apf) 


kinetic energy 


E,(t) = = mut) = =m u? cos? (apf) 
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eo =o 
E,. = max 
fee 


ee 
ae qb 


E. max 


conservation of energy 


lee 
Esmax = Emax m7 Ae 








1 
= —M@, U 
> n =o 


2 Ae 


1 


APPLICATION 


u(x, t) 





The same approach can be used to calculate 
an approximative value of the lowest natural 
frequency of the beam above. 


The method consists in 


- estimating the vibration shape (eigenmode) 


* calculate E(t) and E,(t) 


¢ using E,(t) = E,(t) to get the natural frequency 





Let’stake u(x,t) = ¥ si ri *) sin(ot 


which is the exact eigenmode if M = 0. 
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strain energy 


nO uate 1X 
anes -E vein o r 7) sin(ot 


Ucn} 
al , 2 

Fs = Egpbeam = { >Ellu (x,t)] dx 
0 


= Lely?™ sin? (at) )f sin? ax 
[es 0 = 


[sin Jax = {fi cos{ 27x a 5 


cos(2a) = 1-2sin’(a) 





kinetic energy 


du(x,t) 


u(x,t) = ry 





= = Yaosin| “costar 
ie 


Ex = Exy + Exbeam 
il 2 
Exy = 5 Mah /2,0)] - 


= M Ya" cos? (at) 


m[u(x,t)]? dx 


Nl Re 


Ey beam — 


Co_—_ rr 


L 
1 ny? cos*(at) | sin dx 
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conservation of energy 


L Tt 
Emax =Esmax — o(M+m:) = Elz, 





particular cases 


> | El exact solution since the 
mit eigenmode is exact. 


y 

nm {El EI 
m=0 Oe == 608 = 

2 \ ML? ML? 


In that case, the beam is _ 48El 
considered as a spring with = he 


which gives as 48El _ El 
exact solution @ = pres oars 
ML ML 


M=0 = 


REMARKS 


¢ Only an estimate of the natural frequency 
can be calculated. 


¢ The accuracy of the result depends 
entirely on the shape function which is 


assumed to represent the eigenmode. 


¢ The natural frequency calculated by 
Rayleigh’s method is always greater than the 
exact value. 


SELECTION OF THE SHAPE FUNCTION 


The shape function (eigenmode) must be 
kinematically admissible, i.e. must satisfy the 


displacements boundary conditions at the 
supports. 
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a possibility is to take the deflected shape 


corresponding to the weight of the structure. 
(implemented in some f.e.m. codes) 


example 





CONCLUSION 


The main interest of the Rayleigh’s method 
lies in its ability to provide useful estimation 
of the natural frequency from any reasonable 
assumption of the eigenmode. 





SDF - TIME INTEGRATION METHODS 


p 


P41 
Pj 


Ui lias 


mu+cu+ku= p(t) 


The load p is time discretised. 


P=[ Pp Py Pr --- Pit P; Pist -- Pr | 


The purpose is to calculate u (and y), i 
if required) at the discrete time instants. 


= UO Uh ee Ue Or eaten 
u = [ Up U; U> sas Uj_4 U; Ui) mele | 


p= ad eee etn eee 
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The response at time step /+1 is calculated 
from the equation of motion, a difference 


expression, and known responses at one or 
more preceding time steps. 





Equations of motion at time /and /+1. 
mu;+cu;+ku;=p; (1) 
Mui +0 Uj +K Uj = Pig (2) 
An explicit algorithm uses a difference 
expression of the general form 
Ujen PU Uy Up Ua, Ujey 


which is combined with equation (1) 


An implicit algorithm uses a difference 
expression of the general form 
Uist = F (Uj, Oia Uj, Uj, Uj, Ujy --) 


which is combined with equation (2) 


CENTRAL DIFFERENCE METHOD 


(explicit algorithm) 


difference expressions : 





FPO es elem , 
i : oe assumption 


u=at’+bt+c 





The displacement curve is assumed to be a 
parabola between 3 consecutive points. 
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Introducing the two difference expressions in 
the equilibrium equation at time / 


mu,+cu;+ku; = p; 


gives the unknown u;,, as 


—__ -+ —— 
(At)? 2At 


2am 





initialisation 


The given initial conditions are Uy and Up. 


But Uy and u_,; arerequired to determine u;,. 


The difference expressions for /= 0 are 


ie e 
oy.) ; (At)? 


which after substitution gives 


U_; =U, —Ar uy + 





Up 


(At)? 
2 


The equation of motion at time t= 0 
gives the initial acceleration as 


Py — C Uy —K Up 


Hee 
m 


Ute i 2 y= 2g Us, 
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Algorithm — central difference method 


1 initial calculations 


: = C ly - A 
yy — aS es Uy = Up — Ar dy + Up 


ae ae a tell enc ae bak-2™ 
(At)? 2At (At)? 2At (At)? 


2 calculations for step / 


D; = Pj -au;_,— DU; Ujz1 = Bj/k 


ed ee Pe oie 
if required uy = “Le ee alan ET ie Biel 


3 repetition for the next time step 


replace ;/by /+ 1 and repeat step 2 





AVERAGE ACCELERATION METHOD 


(implicit algorithm) 


difference expressions : 


' Ni ae 
Uj,, =U; + a + Ui) (3) 


(At)? 
A 


Up =U; +; At+ (U;+Uj41) (4) 





The acceleration is assumed constant 
between ¢, and ¢,,. itis given by 


ii(c) = (di +i.)  O<t<At 


The velocity is then linear and the displacement 
quadratic, which gives equations (3) and (4). 
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| Urs 
| Ui 
(3) 


Acceleration . lee e 
(constant) u(t) = 5 ui of Ui) 


i 





Velocity 
(linear) 






eae 
U(r) = Uj + (Uj + dss) 


U; # 


er ne 
u(t) =u, + U; r+ (Uj + 41) 


Displacement 
(Quadratic) 





The system constituted by equations (2) (3) (4) 
has 3 unknowns, ij,, ,Uj,, and u;,, - 


The solution of the system can be calculated as 


2 
Uj) = ar itt — Uj) — Uj 


- 1 . 
Uj.) = ry (Piet = CURA k Ui41) 


The given initial conditions are U)y and Up. 


Py — C Uy — K Up 


liy canbe calculated by Up = a 
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Algorithm — average acceleration method 


1 initial calculations 


ij _ H~CUy—-kKUy pene ok 4m 
: m At (At) 


2 calculations for step / 


“a, 


. il 
Ope = 7 Pit! SC Ue Hele) 


3 repetition for the next time step 


replace /by /+ 1 and repeat step 2 
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STABILITY 


The central difference method will blow up, 
giving meaningless results, if the time step 
chosen is not short enough. The specific 
requirements for stability is : 


en ea 


x 
—— 





The average acceleration method is 
unconditionally stable, which means that the 
procedure leads to bounded solutions s 
regardless of the time step length. 





example 
mu+cu+ku=sin(2z ft) 


@,=8" k=1 &=005 


central difference method: At<0.0796 


ee 
° 





HOW TO CHOOSE At ? 


At must be small enough to get a good 
accuracy, and long enough to be 
computationally efficient. 


For SDF systems, stability criteria are not 
restrictive since At must be considerably 
smaller than the stability limit to ensure 
adequate accuracy. 


Stability of the numerical method is however 
important for the analysis of multiple degrees of 
freedom systems where it is often necessary to 
use unconditionally stable methods. 


One useful technique for selecting the time 
step is to solve the problem with a time step 
that seems reasonable, then repeat the 


solution with a smaller time step and compare 
the results, continuing the process until to 
successive solutions are close enough. 
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NONLINEAR RESPONSE 


lf the structure undergoes large displacements 
or plastic deformation, the equation becomes 


mu+cu+f(u) = p(t) 


where the internal force fu) is a nonlinear function. 


The equation for the central difference method 
can be solved explicitly. 


k. ce = p; —f(u;)— fy ce Lee em U; 
iil peal (At)? DAt at (At)? | 


The equation for the average acceleration 
method must be solved by iteration (implicitly) 
since the first term in [] is unknown. 
of 26 ean 
Salat eee: 
OU|i41 At (At) 





Uj, = Pj+i 


2D CONTINUOUS BEAMS 


Free vibration in bending of a beam with 
distributed mass and flexibility is studied. 


DIFFERENTIAL EQUATION 


U (x,t) u(x,t) ? 


4 
eer erie eal 
ae ee ese 


L_* /9%|, ET m(kg/m) 


pdx 


u ( vf) 124) ve Stax) me owrdx 


ts 
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Vertical dynamic equation 


max?! -y- [v1 Max]+pdx 
at ox 


ovu | aV 


ap hay a = etl) 


Rotational static equation (the inertial moment 
associated with the angular acceleration is 
neglected) : 


_oM 
OX 
ee ee eee 
Ox? dx ax? ax" 


Equations (1) and (2) give 





FREE VIBRATION 


The differential equation 


4 2 
pol + moe = 
Ox ot 


can be solved analytically. 


A possible solution is of the form 
u(x,t) = 9(x)- F(t) 
which gives 


EI 6 (x) f(t) + mex) f(t) = 0 


, ae ’ 
with @4)(x) = — f(t) = 
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a, = -——— =a’ (aconsiant) 
(x) EI f(t) 
d(x) — a9(x) = 0 
> 4 
f(t) + w°f(t) = 0 a = = 


The solution is 


f(t) =Ecos(wt) + Fsin(wt) 


@( x) = Asin( ax) + Bcos(ax) + Csinh( ax) + Dcosh(ax) 





E, F are determined by the initial conditions. 
A, B, C, D are determined by the boundary conditions. 


EXAMPLE: SIMPLY SUPPORTED BEAM 


u u(x,t) = o(x)- F(t) 


Bilom 
x 


Boundary conditions 
U y iti 2y 


M=El—5= E101») f(0 


u(x =0,t) = 0 g(0) = 0 

u(x = L,t) = 0 g(L) = 0 
= 

M(x =0,t) = 0 60) = 0 

M(x=L,t) = 0 o(L) =0 


¢( x) = Asin(ax) + Bcos(ax) + Csinh(ax) + Dcosh(ax) 


Structural dynamics 5.3 


g(0) =0 B+D=0 Bo 
=i =, 
oe D0 


oo A sin(aL) + C sinh(aL) =0 
” — Asin(aL) + C sinh(aL) =0 


a = al) 
Csinh(aL)=0 (2) 


(2) + C=0 or a=0 (nomotion) > C=0 


(1) — A=0 (nomotion) or sin(aLl)=0 
Conclusion: the boundary conditions require 


sin@h—0 and B—C¢ = 0 —% 


¢(x) = a? [- Asin( ax) —- Bcos(ax) + Csinh(ax) + Dcosh(ax)] 
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sinfaL)=0 > al=nm => a= 


visualisation of the three first eigenmodes 


@, are the natural circular frequencies = 
¢g, are the natural eigenmodes El 


Mie) a 
QM, = 41 ie 
An infinite number of solutions u, have been found @, (x) = Asin 7) 
a 


u,(x,t) = ¢,(x) - |E, cos(a@,t) + F, sin(a,f) | 





The general solution for free vibration is a 

- eh 
superposition of all the eigenmodes. @; =9T mit fh TO 
u(x,t) = 2 On(X) - [En Cos(@pt) + Fr sin(a,f) | 


=a Beane the init. cond. u(x,t = 0), u(x,t = 0) 


MDF — FREE VIBRATION 


Free vibration of multi degree of freedom 
systems are studied. 


EXAMPLE 1 - TRAIN 


Differential equations 


Newton’s law is applied for each coach 


eee 


mu, =(k)u,+(-2k) uy 
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mu, +2k u, — kK u, =0 
—> 


mu, —k u,+2k u, =0 


The coefficients in parenthesis are determined by 
considering first u; alone and then u, alone. 


U; > 0 Up = 0 
ku, ku, ku, 
uy =0 u, > 0 


ku, ku, ku, 


et el 


Mass and Stiffness matrices 


The differential equations can be written as 


mitt [oe axl teat 7 f 
/ \ 


mass matrix stiffness matrix 


[m ] [k] 


The stiffness matrix can be determined by the 
classical way. As example, the components of 
first column of [kK] are the external forces that 
must be applied to the masses in order to 
keep the system in equilibrium when a unit 
displacement u,=0 is applied. 
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2k 
first column of [k] = 


u, =1 u, =0 
2k k 


MN > A] <A 


—k 
second column of [k] . 


u, =0 u, =1 
k 2k 


pA] =< AS) 


Eigenvalue problem 


F mot * [ee onl = {th 


The solution of this system is of the form 


a = 12 sina (3 Ont 
U> Py Py 


Since ¢, , @ , a, are constants, derivation gives 


By introducing the two expressions above 
in the system, it is obtained 
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=O, [rm] |2} sinc + (1) 2 sine = (3 


+ ([ie}- fm} |*'} sino,n = 14 


2 


This equation is valid for every t, which implies 


eigenvalue 
problem 


[1k] [m]) |$ 





This system has one trivial solution (¢, = @ = 0) 
which corresponds to equilibrium (no motion). 


Other solutions can be found if the following 
condition is respected: 


det ([k]- 3 [m]) = 0 (2) 


(2) > 





Die = co iahs 
en 


+ (2k-a-m)—-k? =0 


Vk am Okami kh) 0 





Two values of a, have been found. For each 
of them, 6, 5, can be determined by 
introducing the value of @, in the system (1). 


(1k}-o [m]) 8) = {0} on 


i) 0 
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(am =k) 


o+[t fs} 


This system has an infinite number of solutions, 


lal 


one of them is 





Conclusion : Two solutions have been found 


{64} 2 
A ee 
u(t 


SS 


Uj(t)| _ i in(.(3k/m - t) 
i eme 


@, @, ‘natural circular frequencies 


10,} {0} : eigenmodes 





Every linear combination of this two solutions 
is also a solution of the free vibration problem. 
The general solution can be written as 
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Gee) eee Cle 1) 
eee =C, {i sin(a,t+0,)+C, {. ; sin(@,t + >) 


Physical interpretation of the eigenmodes 


Free vibration Is initiated by an initial deflection 
corresponding to eigenmode 1: 


u(O)=A u(0)=A t&(0)=0 4,(0)=0 


With these initial conditions, the response is 


CG) ema U |er 
ie a tl a 


¢ The motion of both solids is harmonic with 
qo, for circular frequency. 


¢ The deflected shape is constant in time 
and corresponds to eigenmode 1. 





EXAMPLE 2 —- THREE STORY FRAME 


m=2kg 
L=200mm 


columns 
b=10mm 


h=1mm 





E=200GPa 


(4 columns by story ) 
48El bh* 
—— ——— 
b 2 


me of © 
[kJ=|=k 2k =k 
(ames 


[m] = 


oo 3 
eS 2S 
3.00 


( [k] is symmetric ) 


u, =1 — k 


< k 


<k 
u, =1 — 2k 


<k 


u, =1 — 2k 
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first column 
k 
—k 
0 


second column 


—k 
2k 
—k 


third column 
0 
—k 
2k 


Eigenvalue problem 
([k]-«? [m]) {0} = {0} 
This is solved numerically by MATLAB. 


M =[m| 


[F, E] = eig(K,M) 
K =|k| 


solves the eigenvalue problem [k]{o} = A[m]{o} 
> @=A 


The MATLAB solution is given as 


oO, Oe 
F=| 0 O; 0 F =| 4017490 7493 
EO 
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The eigenmodes are determined with an 
arbitrary multiplicative constant. By convention in 
dynamics, the eigenmodes as given on the form 


Pi 
0, = \% 
P31 


The eigenmodes in MATLAB are calculated such 


D1 
Oo = \P1 
P31 


with 0: +051 +9, =] 


In order to obtain the dynamic convention, the 
following scale command must be used. 





F = [F(:,1)/F(11) F(:2)/ F(12) F(:,3)/ F(1,3) | 


Results given by MATLAB 


@, = 9.95 rad/s @, = 27.9 rad/s 
f, =158 Hz f, = 4.44 Hz 
1.000 1.000 
o, = 40.802 b, =4—0.555 
0.445 — 1247 


=() 355 


— 1.247 





Physical interpretation of the eigenmodes 
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@; = 40.3 rad/s 


1.000 
3 = 4-2.247 
1.802 


aA) 


1.802 





If free vibration is initiated by initial displacements corresponding to the eigenmode n, the 


vibration of each story will be harmonic with frequency f, and the structure will vibrate with a 





constant deflected shape corresponding to the eigenmode n 
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Comparison between discrete and distributed systems 


the eigenmodes 
discrete systems are defined with 
a multiplicative 


constant 
eigenmode = vector 


number of natural 
frequencies = number P| 
of degrees of freedom o, = +, 


P31 ieee! 


convention 





distributed systems eat eee eigenmode = function convention 


natural frequencies 


A=1 
Pre] ee Asin 


MDF — MODE ANALYSIS 


The mode superposition approach is used in 
order to calculate the response of a MDF 
system to an applied load. 


First harmonic loads and then arbitrary loads 


are considered. Finally damping is introduced. 


HARMONIC LOAD - TRAIN 


The first coach is excited by a harmonic load 


U; U> 


— 
Py sin(at) 


Ce ed 


= 
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Summary of the eigenvalue problem 


eigenvalue problem ([k]- @? [m] ] {o} = {0} 


a, =k/m {oJ= {i Q, = ./3k/m {o.}= {|| 
The matrix [o] is defined as 


l-[feHotl=|) 


Differential equations 


lot Ma a] Se 


mu, + 2ku, — ku, = p, sin(at) 
Vi aU 2K Ua) 


The two equations are coupled and it is hence 
not possible to solve directly the system. 


Modal coordinates 


The idea is to change coordinates (unknowns) 
and calculate first g,(f) and q,(t) which are 
defined by 
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by derivation, it is obtained 
U,(t) q(t) 
to } = ol (2 
u,(t) q(t) 
Introducing these two expressions in the 


differential equations and multiplying on the 
left side by [o]' gives 


fol" [m] (ol) +o el} = [ol {p}sin(oo 


which can be rewritten as 


[M] | + [k] a = {P} sin(at 


oD) 2 


[M] = [o]' [m] [9] 


{P} = fol" {p} 
(K] = fol” [k] [ol : 





After calculations it is obtained 
2m 0 2k O 
[M] = [kK] = {p}= 4" 
O27 0 6k Do 


age 0) a i" 2h OG; = ps en 
QO 2M} (q2 0 6k} \% Po 

2mq, + 2kq, = Py sin(at) 

2Mq, + 6K q> = Pp sin(at) 


Since Both [M] and [K] are diagonal matrices, 
the two equations are uncoupled and can be 
solved separately. 


Each equation represents a SDF system. 


The natural circular frequencies for the first and 
second equations are @ and @, which have 
been calculated by the eigenvalue problem. 
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Each equation are solved by considering 
only the steady state response 


a Po /2k nee ie 
q(t) ae sin(at) =a = k/m 
qo(t) = aie sin(@at) @, =/3k/m 

1—(@/a,) 


Uy; and Us can now be calculated with 
ui(t)| — 
ne e 


Ut) (2h — | sin(at) 


l-(o/a)y 1-(o/a,/y 


Po 2k _ __ Po/6k | sia 
1-(v/a) — 1-(@/a,) 





Conclusion 





_ | Po/2k Po/6k | .. 
ae fem ‘ errs] a 


yes) 2s CK 
Sl (eres dares 





After some time (steady state response), the 
structure vibrates with the same frequency as 
the force. 


Resonance 


the load frequency @ is close to the first 
natural frequency of the structure @, 


al ae 
oO, & 1 (a) P sin(at) 
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¢ the structure vibrates with a deflected 
shape corresponding to eigenmode 1. 


¢ the amplitudes of the vibrations become 
infinite (undamped case) 


lf a MDF structure is excited by a harmonic 
force whose the frequency is one of the 
natural frequencies of the structure , then 
after some while (steady state response) 


¢« the structure vibrates with the same 
frequency as the applied force 


¢ the deflected shape is the eigenmode 
associated to the natural frequency 


¢ the amplitudes of the vibrations are large 


This phenomenon is called resonance 
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HARMONIC LOAD — 3 STORY FRAME 
Applied harmonic force 








A harmonic force is applied to the second story. 0 
1‘ sin(2xt) = {p} sin(2z%) 
m 
as — wae 0 
Im] =/0 m 0 
Lsin(2zt) Oo Oo m nee 
a > vu, Mode superposition 
oe u(t) qy (t) 
ee > LU, [k]=|—k 2k —-k u(t) = lo] 4 q(t) 
0 =k 2k U3(t) q3(t) 
Uy Uy 
Eigenvalue analysis [m] .a$ + [k] u,+ = {p} sin(2zt) 
U U3 
| (9.951 1.000 1.000 1.000 ; 
wo, = 27.88$  [o] = | 0.802 -0.555 —2.247 a1 q | 
40.29 0.445 -1.247 1.802 > IM] sq) + [K] jq/- = {P} sin(2nt) 


ws i 
Q3 93 


3.682 0 0 
[M] = [o]' [m][o] =| 0 5.726 0 
0) 0) 18.59 


364.7 0 0 
[K] = [o]' [k]lo] =] 0 4452 0 
0 0 30184 


0.802 


{P} = [o]' {p} = }-0.555 
=—2. 747 


The uncoupled system is 


3.682 G, + 3647 q, = 0.802 sin(2z?%) 
5.726 G, + 4452 q, = -—0.555 sin(2z?%) 
18.59 G, + 30184 q,; = —2.247 sin(2nt) 
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Each SDF system is solved by only considering 
the steady state response, which directly gives 


gy (t) = 280273087 sin(2nt) = 3656 10-3 sin(2nt) 
1—(27/9.95) 
,(f) = 0? sin(2mt) = — 1313104 sin(2nt) 
1— (27/279) 
g,(t) = 4100184 inant) = — 7.629 10° sin(2nt) 


~ 1-(27/40.3) 


Transformation to real coordinates 


u, (t) q,(t) 3.656 10° 
u,(t)+ = [o]4q,(t)$ = [o]4 -1313 10-7} sin(2xt) 
u;(t) q;(t) SoZ Ome 





u, (t) 3.449 10° 
u,(t)$ = 43177 10°} sin(2xt) 


u,(t) 1.653 10° 





ARBITRARY LOADING — TRAIN 
p (N) 


p(t) is applied to 


the first solid bs 


0 1 
mf Mah = Co 


Transformation to modal coordinates 


cw ‘ahs ik] 4) te 
> Qo 


2mq + 2kq = ptt) 

= ee 

2M q, + 6K Qq = Plt) 

Each SDF equation is solved by using a time 


integration method (lesson 4). For this particular 
case, the exact solution can be found (lesson 3). 
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1 


‘ N Gq =leeeee | 
— g, =[eecce| 
0.4 

0.2 4 

0 
-0.2 
-0.4 
-0.6 
-0.8 

A t (Ss) 


0 0.2 04 0.6 0.8 1 fee 1.4 1.6 1.8 2; 


Transformation to real coordinates 
oy ae q(t) 
(oo - i (a 


tor a eee e ee . eee 
U, (ft) eee ae Q), eooc5e 


must be performed 
for each time step 
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SYSTEMS WITH DAMPING 
Solution : 
The equation system of MDF 


ae ee 1) Natural frequencies and modes are 
systems with viscous damping is 


calculated by neglected damping. This is 
correct if the damping is small. 
[m]{d}+[e]{uj+[k]{uy = {ple 
2) Transformation to normal coordinates : a 
The transformation to modal coordinates damping coefficient €, is introduced for each 
can be performed as before, which gives mode. & is determined by experiment. This 
gives n uncoupled equations on the form. 


tuy=loltat [M]tqy+[C]iat+[K] ia} = {P(A} 
with [C] = [9]' [elo] 


Ch 
2,/K,M, 


Mp Gn + Ca Qn + Kn Gp = F(t) cn= 


3) Each SDF equation is solved by time 
integration method 


—> {q(t)} 


Problems : 


- how to define [o] ? 


4) Transformation to real coordinates 


tuy=loltay 


¢ how to obtain by experiments [C] or [c] ? 


¢ useful transformation only if [C] is diagonal 





DAMPING MATRIX 


For certain problems, it is better to solve directly 
the coupled differential equations 


[m]{uj+[e]{ut+[k] {uy = {p(t} 


instead of using the modal transformation (see next 


section). The damping matrix [c] is then needed. 


[c] is calculated from the modal damping ratios €, 
which are determined by experiments. 


Method one 


[c] is calculated from the modal matrix [C] 


__ G, 2 
C3 


[c] = [ol [e] [6] > [el = (fo) [ello] 


This method supposes that all the damping ratios 
€, are known, which is usually not the case. 
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Method two : Rayleigh damping 
The damping matrix is taken as 
[ce] = ay [m] + a [k] 


It can easily be shown that the modal 
coordinates transformation leads to a diagonal 
modal damping matrix and that 


Hence, [c] is calculated using only two modal 
damping coefficients. 


The two modes with specified €, should be 


chosen to ensure reasonable values for the 
other damping ratios. In practice, the lowest 
modes and the third or fourth lowest ones 
are used to determine a and a,. 





EFFECT OF THE DAMPING ON 
HARMONIC VIBRATIONS — TRAIN 


The first solid is excited by a force p, sin(at) 
Transformation to modal coordinates 
2mMqG,+ C, Gg, + 2kKq = Pp sin(at) 
2Mq,+ Cy Q + 6KQ = pp sin(at) 


The steady state response for each equation is 
on the form 


= D,sin(at-@) 
q(t) = D, sin(@t—@,) 


Py /2k 


= . 
h-(@/@)? |} + bg (o/@)]? 





ns aa 
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Transformation to real coordinates : after 
some work, the response can be written as 


ee : oa i ote 
U,(t) q(t) G, sin(at — a) 


Conclusion : 


¢ Both solids vibrate with the same frequency 
as the applied force, but with a _ certain 
difference of phase. It means that u,(f) and 


u,(t) don’t reach their maximal values at the 
same time. 


¢ The resonance properties enounced without 
damping are still valid. 
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MDF - TIME INTEGRATION METHODS 


CENTRAL DIFFERENCE METHOD 
The system of differential equations can also be (explicit algorithm) 
solved by time stepping methods. These methods 


are similar to the ones for SDF systems (lesson 


This method is based on two difference 
4) by replacing scalars by matrices and vectors. 


expressions and the equilibrium at time /. 


ral tureerel ee cent piel? fu} ee fg 


2At Ty 
Initial conditions  {u},  {u}, 


[m] {uj at [ce] {uj a [k]{u}; = {Pp} 








combination of these 3 eqs. gives - Ty lt uj;., = {Pp}; - -| a El tu Uji_1 - i a 


2 2At 





(At)” 
using the following approximation = lel {tu _-{uh.,) = ale |({u}_, -{u}) 
it in 2 : Im] _ [c] 2|m] | [c] 
tis obtained pbs = (0h | — Lf, |e) AF uy 


Then, if [m] is diagonal, the solution of simultaneous equations is not required. 


AVERAGE ACCELERATION METHOD 
(implicit algorithm) 


For this method, a set of uncoupled equations 
cannot be obtained by simplification since the 
equations are on the form 





ic) 2 A fu = {p} 


At (At) 


And if [m] can often be taken as diagonal, this 
is not possible for the stiffness matrix [k] 


CONCLUDING REMARKS 


The modal superposition requires that the 
system Is linear, i.e. small elastic deformations. 


By comparison to time integration methods, the 


main interest of the modal analysis is that in 
most of structural dynamics problems, only a 
few modal contributions need to be solved. 
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Example 
[m]{u}+[e]{u}+[k]{u} = {p(t)} 10000 d.o-f. 


— 10000 modal equations 


Mn Gn + Cn Wn + Kn Qn = F(t) 


Only the 10 first equations, corresponding to the 
10 lowest natural frequencies need to be 
solved. 11 .--. 910000 can be neglected and don’t 
need to be calculated. 


U, (t) Q, (t) 
uot) | — tq | aol) 
U, ; (t) = lol 0 
Uroo0o() 0 | 


{u(t)} = {,}.q;(t) 7 {> }q>(t) ve + {10} Holt) 
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Implicit time integration methods are suited when many modes would be needed in the modal method. 
These methods are often unconditionally stable and the time step At is only limited by accuracy. 


Explicit time integration methods are often used in combination with diagonal mass matrices which 
allows a small cost per time step. These methods are often conditionally stable and the time step must 
be very small so that stability is ensured for all the modes: even if the response in the higher modes is 
insignificant, it will blow up if the stability requirements are not satisfied relative to these modes. 


Contrary to implicit methods, explicit methods do not require an iteration process for nonlinear problems. 





INTERPRETATION OF THE MODAL SUPERPOSITION 


u, (t) Q(t) Pi Pr) P31 

Un(t)> = [ol 4 ao(t)p =4A2F O(t) + 426 Q(t) + 432° Q5(t) 

u3(t) q3(t) —3 Pr3 P33 
eigenmode 1 eigenmode 2 eigenmode 3 


= q(t) + (0) = G@.(1) 


EARTHQUAKE ANALYSIS 


An introduction to earthquake analysis is given. 


Analysis by response spectra is presented. 


SDF SYSTEMS 
ASDF system is subjected to a ground motion 


u,(t). The deformation response u(t) is to be 
calculated. 


m(Ug +U)+cu+ku=0 


U+2éEo, U+ a, U=—i,(t) 
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The ground acceleration can be registered 
using accelerographs, see the example below. 


A ht} aM Nil) 


0.1 | 
0 HW i pa HAH Ht TW he ' 


| V7 kth A 
ven ‘ail 


Acceleration 
Acceleration of gravity 


0 5 10 15 20 25 
Time, sec 


FIGURE 24-15 
Accelerogram from El Centro earthquake, May 18, 1940 (NS component). 


EQUIVALENT STATIC FORCE 


| +| u(t) f(t) f(t) = k u(t) 
=m A(t) 


2 = ma; u(t) 
: 
| 


A(t) = w-u(t) # u(t) 


Pseudo accelerati on 


f,(t) is the force which must be applied statically 
in order to create a displacement u(t). 
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REPONSE SPECTRA 


A response spectrum is a plot of maximum response (e.g. displacement, velocity, acceleration) 
of SDF systems to a given ground acceleration versus systems parameters (T,,, ¢). 


A response spectrum is calculated numerically using time integration methods for many values 
of parameters (T,,, ¢). 





Example : Deformation response spectrum for El Centro earthquake 





10 
T, = 0.5 sec 0 
C=2% 
» 8 -10 
3 10 
‘3 & 
T, = 1 sec & 
(=2%| & ° s 
ee 
a A -10 5.97 in. Q 
10 
T, = 2 sec 0 
Ca2% 
7 -10 7.47 in. 
0 10 20 30 





Time, sec 


Deformation, pseudo-velocity and pseudo- 
acceleration response spectra can be defined 
and ploted on the same graphs 


Peak Deformation D = max|u(t) 


Peak Pseudo- velocity Vo 





Peak Pseudo-—acceleration A=a- D 


q, : natural circular frequency 
of the SDF system. 


Hovity V, in./sec 


Pseudo-velotity V, in./sec 
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COMBINED D-V-A SPECTRUM 





EXAMPLE 


A water tank is subjected to the El Centro 
earthquake. Calculate the maximum bending 
moment during the earthquake. 


m = 10000 kg 
7 k = 98.7 kN/m 
E=2% 


oy =| =3.14 rad/s > Te ee 
m 


10m 


L= 


Qn 


D=7.47-25.4=190mm 


Spectrum > 
A=0.191-9.81=1.87 ms 


(obs: A=; D) 
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V, in./sec 


0.2 
0.02 0.05 0.1 0.2 





100 


50 
A AN Aronm fn. 


f =k-D=18.7kN 


Moray =187kKNm 


When the equivalent static force has been 
determined, the internal forces and stresses 
can be determined using statics. 
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RESPONSE SPECTRUM CHARCTERISTICS T, = 2m mlk 


General characteristics can be derived from 
the analysis of response spectra. 


T, < 0.03 s : rigid system 
no deformation 





u(t)=0 >D=0 
— 
T,>15 8s: flexible system 
no total displacement 
u(t) = u(t) > D= Ugo 
= 


The spectrum can be divided in 3 period ranges : 


t= 05s : acceleration sensitive region 


0.5<7,<3s : velocity sensitive region 


f= 3s : displaceme nt sensitive region 
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ELASTIC DESIGN SPECTRUM 


EXAMPLE 
Problem: how to ensure that a structure 


will resist future earthquakes. 


Elastic design g Aylty., d 
spectrum 





The elastic design spectrum is 
obtained from ground motions data 
recorded during past earthquakes at 
the site or in regions with near-similar 
conditions 


Peak ground acceleration, 
velocity, and displacement 


Pseudo-velocity (log scale) 





10 sec 33 sec 
1/10 Hz 1/33 Hz 


Natural vibration period (log scale) 


Figure 6.9.3 Construction of elastic design spectrum. 
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MDF SHEAR BUILDING 


m 0 0 (itu, ti u) (0 
O m, 0 |i, +g; + [e],ty> + [k]{uy¢ = 40 
UO ia Uy + Ug U3 U3 0 











m 0 O U; U; Uy; Mm 

O my 0 [stp + [ce]; a> + [k]4u ¢ =-4 mM Olt) 

0 0 m| |U; U3 U3 Mm, 

Eigenvalue analysis dh bt. 03 
Deets > = 2 0» “|e 03 = es lo] = & Py es 


P31 9) P33 P31 P32 933 
Modal transformation 


i C, q A +|™ 
Cos Cy ou =—[o]' m Ug(t) {P}=|[o}' 5m, 
93 C3 | (93 - Ms 
Mn Qn + Cn Qn t+ Kn Wn = — Fy Ug (t) f= 23 
r,=— 
. ; . M 
Gn + 26n Wy Gn + Op Gn=—Tln Ug(t) z 


M, Ky 


M, 
M; 


Ky 
k; 
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RESPONSE HISTORY ANALYSIS _—_( i/,(t) known) 


Solve numerically Gn + 2€ pn In + Op Gp =p lig(t) n=1,23 > q(t) n=1,2,3 
Summation of the u, (t) Gy (t) —1 P~i2 —13 
modal contributions Uy(t) > = [o]s ar(t) > = 4dr p HH +4 O22 f H(t) +5 3 7 Blt) 
u3(t) q3(t) P31 P39 P33 
Equivalent static forces 
Observation : The f(t) > — u(t) 
f(t) u(t) g,(t) contribution of the 
f(t) =[k]Jun(t) =[k] fo] a3(0 Tne ymoee> sea f,(t) > —> up(t) 
f,(t) u,(t) g(t) often be neglected. : ‘ 
it f(t) > —> U;(t) 
f i = [k] a0 q(t) + [k] a y(t) + [k] es q3(t) 
F(t) P31 P39 P33 
V, (t) 
The internal forces and stresses can Example : base shear force 


now be determined using statics. 
static V,(t)=f(f)+A(t)+6(t) Vemax =max|V,(t) 


RESPONSE SPECTRUM ANALYSIS 
( Ug(t) unknown) 


The design response spectrum is used to 
estimate the behaviour of the shear building. 


The modal transformation gave (page 7.7) 


n= 123 


Gn + 26n On Qn + Dp Qn=—T pn Ug(t) (1) 





which can be compared to the equation 
for a SDF (page 7.1) 


+2, 0, U+@7 u= —Ug(t) (2) 


D,, , the maximum displacement for equation (2) 
is found on the design response spectrum. 


The maximum displacement for equation (1) is 


eS I, Dy 
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Summation of the modal contributions 


Ui max ry D, 
but Uy max = lo] I, D, (3) 
is 


U3 max 


because D,;, D2, D3 do not occur at the same time. 


(3) over estimates the maximal real displacements. 
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The summation of modal contributions is performed in a different way. 


The idea is to calculate first maximum displacements and equivalent static forces for each mode. 





hm) > — Uym1) h(m2) > — U(m2) 

Fm am —_ U3(m1) Bom 2)—> — U3(m 2) 
Ui(m1) ,D, Pi Ui(m2) 0 V2 
Umi = lol} 0 + = 44/7 MD, Usim2) ¢ = 114 Dy = 402 ¢ Dy 
U3(m1) 0 P31 U3(m 2) 0 P39 
f(mi) lyD, Yi him2) 0 M2 
him) ¢ = [k]lo]; 0 }=[k]i a, ) 74D, fm) p= LK] lo] 5 2D, > = [k]; 5 - FD, 
Fm1) 0 P31 Fam) 0 P39 


Same calculations for the third mode 


Once a value of the required response has 
been calculated for each mode separately, the 


superposition of the different modes is 
performed using an approximative rule. The 
most common one is the SRSS approximation. 





Maximum displacements 








= ‘ U5(m1) a U5 aval U5 3) 


, U5. one U5 (ry Dhak U5 (my 3) 


Uy max 





U ‘ Ui(mi) + Ui(m2y + Lima) 
| lmax | 


U3 max 


Maximum shear force in the highest floor 


af 
Yi 
<j} ————_. 





Vimax = fa (is Aim ae fim 3) 
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Maximum base shear force 
Maximum base shear force due to the mode 1 
Vomt) = Amt) + foamy + fm1) 


Maximum base shear force due to the mode 2 


Vim2) = fy (m2) + f(m2) + m2) 
Maximum base shear force due to the mode 3 


Vim 2) = 4m) + fom) + (m3) 


Total maximum base shear force 





2 2 2 
Vomax = Af Vicm1) + Voim2) + Vom) 


Observation : The contribution of the 


higher modes can often be neglected. 
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FINITE ELEMENT METHOD 


The generalised load vector and stiffness 
The equation system obtained by f.e.m. is matrix have been studied in the f.e.m. 
course and are identical in static and 
dynamic analyses. 


[m]{u}+[e]{uj+[k]{uy = ip(o)s 





The damping matrix is not obtained by finite 
element discretisation. The different ways of 
introducing and dealing with damping have 
been studied in lesson 6. 


with 


mass matrix The resolution of the equation system, by 


mode superposition or by direct time 
integration methods has been studied in 
lessons 4 and 6. 


damping matrix 


stiffness matrix 


The only remaining work is to present the 
derivation of the element mass matrix. The 
assembly of the element mass matrices to 
the structure mass matrix is performed in 
the same way as for the stiffness matrix. 


generalised load vector 


displaceme nt vector 





CONSISTENT ELEMENT MASS MATRIX 


The equation of motion can be derived by 
equating the work done by the _ externally 
applied loads (external work) with the work 
absorbed by inertial and dissipative forces 
(internal work) for any virtual displacement (that 
is, for any imagined small motions that satisfies 
compatibility and essential boundary conditions. 


W; =| ({8e}" {o}+{8u}" p {u}) av 


i=l 





We = [{du}" {F}dv + |{su} {@}ds+ > {su} {ft} 


where {F} and {®} are the prescribed body 
forces and surface tractions, {f}; and {du}, 
represent prescribed concentrated loads and 
their corresponding virtual displacements at a 
total of n points, o is the mass density. 
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{du} and {de} represent virtual displacements 
and their corresponding strains. 


Finite element discretisation provides 


{u}=[N]{d} {u}=[N]{d} {e}=[B]{d} 





Shape functions [N] are functions of space 
while nodal d.o.f. {d} are function of time. 


A linear elastic material is assumed 


The internal and external virtual works can 
then be rewritten as 


Ww, = {Sd}! | ((@1 {elte1{a1~ (NT n)fa}) ov 


V 


W, = {8d}" FINI {F}av+ J [N| {@}ds+ > Nite, 


iI 


The last equations in the preceding page are 


Ww; ={8d}" | ([e)"[e][B]{d}+[NI" o[N]{a}) av 


V 


i=l 


W, = {Sd}' | IN] {F}dv+| [N]'{@}ds+ > [NJ {f} 


W,= W, is true for arbitrary {5d}, which gives 


fe IN] [N] dv {d}+ f [B][E][B]dv{d}= {p} 


{p}=j [NJ {F}av+ [N]'{@}ds+ 3 IN] {f} 


| 


| 
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This can be rewritten as 


[Lm] {u}+[k]{u} = ip(e)s 


with 
[k] = j [B]'[E][B] dv 


[m]= J o [NT [N] dv 


V 


The mass matrix defined in the last equation 
is called consistent element mass matrix. 


The word consistent emphasizes that this 
form has been derived using the same shape 
functions as the element stiffness matrix. 
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BERNOULLI 2D BEAM ELEMENT 
shape functions 


Vi Vo 
6, 8, - 
NU etm ‘ m=pA AN 





1 L 2 o 3 - 
a(x)=1-3[ 7) +2(2) x)= x{ 1-2) 
v(t) 
_ A(t) 
v(x,t) = [N] td} id} = me 
0,(t) ee |: Ny 
x) x) x {|x 
IN] = La(x) 92(x) s(x) 94(x) ] os) =3( 7) -2(#] 9400) = (E-1) 
156 22L 54 -13L 12 6L -12 61 
| —mL| 22k 4P 13k -30 _ Fil 6b 42 -6L 20 
er ee es Lm] = 250 oy an i) A aa ee heen Ome 


Se ie Deere Gil ale Gir Ale 


[m] 


Structural dynamics 8.5 
BERNOULLI 2D BEAM ELEMENT WITH AXIAL DEFORMATIONS 


0, A bar element is superposed to the previous beam element. 
| ean | ‘ 
u(x,t) = Q(X) Uy + P(X) Up 
1 u; iE 2 Uy 
A A 


The shape functions for the bar element are 


p(x) =1- . P(X) =F 
1440 0 0 70 0 0 EA/L 0 0 -EA/L 0 0 
Qin i56n 22 n Ome 4 ior 0 126/28 6EI/P 0 1261/2 6EI/L 
mL | 0 22k 4? 0 13k -30 rk] = 0 66/2 4EIV/L 0 -6EI/2 2EI/L 
420/70 0 0 140 0 0 -EA/L 0 0 EAL 0 
stoi 0 560 


0 
1261/2 —6EI/L 
~6EI/2  4El/L 


0 -126//2-6EI/P 0 
OF se 0 eae 0 6EI/L 2EI/L 0 


LUMPED MASS MATRIX 
For certain applications it is better to use a 
diagonal (or lumped) mass matrix. 


Two obvious advantages are less storage 
space and less processing time, especially in 
case of an explicit time integration scheme. 


PARTICLE MASS LUMPING 
Example : Bernoulli 2D beam element 


Vi Vo 
A, 0; 
> Elim 


2 


m=pA 
: L 





An easy way of obtaining a lumped mass 
matrix is to replace the distributed mass by 
two particles of mass mL/2 at each node. 
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The rotational inertia is defined by considering 
that a uniform slender bar of length L/2 and 
mass mL/2 is attached at each node. 


The associate inertia moment is 


J =(mL/2)(L/2)*/3 


mb 
mL = 1 24 mL «% 
2 Zz 
U; Uy 
1 ss B 2 
24 


i 6 
mL|0 2/12 0 O 
[m] = 
Pai Ome 
@ 90 Oo (ee 
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HRZ LUMPING 


Different methods can be used to transform the consistent element mass matrix and obtain a 
diagonal matrix. One of them is the HRZ lumping. 


The idea is to compute only diagonal terms of the consistent element mass matrix, then scale 
them so as to preserve the total element mass. Since there may be both translational and 
rotational d.o.f., the method for an element of total mass M is 


1. Compute diagonal coefficients m;, of the consistent element mass matrix. 


2. For each coordinate direction in which motion is described by the element d.o.f. 
a. determine a number S by adding the m, associated with translational d.o.f. only. 
b. multiply all coefficients m,; associated with this direction by the ratio W/S. 





Example : Bernoulli 2D beam element 


4. 


[56n 0 any ee: ie Oe BO a eA0) 


; 312 : 
Ooo On 0 S=——m mL|0 7/39 0 O 
[m] = — 420 Im] = = 
250 OntoG nO vane 7m cee nO) 
0 0 0 4e as 0 0 0 27/39 
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EXAMPLE 


Mass and stiffness matrices for the structure 
The natural frequencies and eigenmodes 


of a simply supported beam are studied. Le=L/2 
The beam is discretised in 2 elements. . 
io eer =) rn) 0 
laa Gk ie ts 0 
EIALm ik] _ il -12 - GLa 12+12 = Gla OLn ! —12 6L, 
fs A Hoey ce erate een eee alee tes 
en esc Deer 
m= p A ! i: i 
0 0 Cen ie enema 
A, a) 03 
{ — gE 54 -13lg | 0 0 
2 | 
;, t ees 13L, 22 sO 0 
V> V. Roe eees eee - SSS = Se 
" Ty; bm) = Mee 54 13lg «1564156 |- 2254221, 54 — 13Lg 
= 2 2 ae 2 
ee ee 420 |-13lg -314 -22lgt2le | 41G+41G | Ble -315 
oO | SA ice 1s = ks 
0 0 TS cele =) Leia 
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Boundary conditions 


(P,(t) , 23(t) unknowns 





Equation system 


Ale Eee ea On 6) Ai =o 2 0 |e te 
rr) ele ee OTe 00 faeces tC) 
420 =i ee PS eee, ie Di aU DeClE meee: |G) ~ Jo 

es ee see ema): OrOle 2 are.) eae 
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Solution of the eigenvalue problem (exact values in parenthesis) 


Ej EI} El E} 
@, = 9.91 ,|—— W, = 43.8 ,| — @3 =110.,|—— @, = 201, |—— 
ml ; ml : ml : ml! 
(9.87) (39.5) (88.9) (158) 
] I ] ] 
oe 032L eens ae —0.05 L ee 
| 20 me ea lo ae 
—] ] —1 ] 
eigenmode 1 eigenmode 2 eigenmode 3 
0; = 0 





EXPERIMENTAL DYNAMICS 


The purpose of experimental dynamics is to 
determine by measurements the dynamics 
characteristics of a structure. This is performed 
by measuring the response of the structure to a 
known excitation. 


known known 
excitation response 
unknown dynamics 
——————_> ————_—> 


characteristics 





Unknown dynamics characteristics 


SDF : natural frequency  f, 
damping coefficient ¢ 


MDF : natural frequencies f, f, f, ... 
eigenmodes 0, O O,... 
modes damping Gece mcane. 


Structural dynamics 9.1 


Different types of excitations can be used 












Free vibration : the structure is disturbed from 
equilibrium and then vibrates 
without any applied forces. 


Harmonic force : a sinusoidal force is applied 
through a special device. 


Pulse force : a hammer is used to apply a 
short pulse force. 


Random force : random noise as e.g. traffic or 
wind loads are applied. 


For practical reasons, accelerations are 
easier to measure than displacements 


SDF SYSTEMS 





Objective : determine Cc 


experimentally f,andé k py 


Free vibration test 


u(t) = Ce >”! sin(a,t + 6) 


+> i(t)=De°™ sin(a,t +a) 


After derivation a similar expression for the 
acceleration is obtained and consequently 
the results derived in Lesson 1 for the 
displacements can be used. 
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Harmonic force test 


The structure is excited by a harmonic force 
whose frequency is slowly increased step by 
step in order to reach the steady state 7 
response at each increment. The amplitudes 

of accelerations are measured at each step. 


Results from Lesson 2 : 
p(t)=p,sin(at)  u = |u| sin(at+ 9) Ae 


PK 
Wh-(olo,)27 + [2E(@/a@,)] 





ul = 





In reality, the maximal value is not known and 


By derivation, it is obtained : 
a curve fitting is done. 


ui = |i] sin(at +) i] = oO” ul lf the damping is law, the same method can be 
used with lu instead of lu. 


Ambient vibration test 


The structure is excited by a small random force, 
which means that the force is not known. 


This method is often used for bridges where wind 
and traffic loads are used as random forces. 


uae Ceo sin(@,t+a@) + Particular solution 


A ~ 


initial conditions Ug Up random load 


u a 
a es 
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The response (acceleration) is the sum of the 
free vibration solution and the particular solution 
which depends on the random load. 


The idea of this approach is to remove the 
particular solution and the initial condition yj, by 
taking the average of many subrecords of same 
length At and same jj, . The result is a free 
vibration solution which can be then used. 


The trigger value (j, and the length At must be 
carefully chosen . 


At 


_h kT nen 
FIA AL EVV) PA 
SO ee 


MDF SYSTEMS > U 

: i : => Up 
Objective : determine 

experimentally f, ft Se 


; 2 3 and ¢; op o3 


Harmonic force test 


The structure is excited by a harmonic force 
whose frequency is slowly increased step by 
step in order to reach the steady state 
response at each increment. The amplitudes 
of accelerations are measured at each step. 


Results from Lesson 6 : 
U, (t) O1 Ar A3) A, 


Uy(t) > = | by br th3 |4 Ao 7 SiN(@t + 8) 
U;(t) ~31 Pry 33 | (A; 
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On =) 
Uy (t) 1 
A, large 
am aaa U, (t) = D>) A, sin(at + 0) 
A, , A3 negligible ; 
u(t) Pr} 


The structure vibrates with a deflected shape 
corresponding to eigenmode 1 and the 
amplitudes of vibrations are large : resonance. 





The sign of ¢, is determined by the phase 


which is also measured. 





ii, 1.00 
Uy p =C.-0.55> sin(@t + 30°) 
alee) 


Explication : the 
solution at the first 
resonance is e.g. | ti, 


1.00 sin(at + 30°) 
sullen = C { 0.55 sin(at + 30° + 180° 
rewritten as }2/ ~ ee) 

Us 1.25 sin(@t + 30° +180") 


The measured phase is either 30° or 210°, 
which gives the sign of ¢; . 


1, $2, $3 are determined by band width method 


for each resonance. 
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Remark : this method assumes that only one 
mode is acting at each resonance. For that, the 
natural frequencies must be well separated: 


@ =a; — A; large and other A; negligible 


Aj 





In that case, the influence of the 2" eigenmode 
can not be neglected at the 3'¢ resonance. 
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Free vibration test (damping neglected) 


Uy, (t) 1 02 3 
Results from Lesson 5 : Un(t)> = 41; SiN(at+,) + +» - Sin(@,t+4,) + 437 SIN(@;f + A3) 
u(t) Ps} P30 P33 


U;, (t) = Di sin(a,t + 6) oh D> sin(@,t + 05) ie 23 Sin(@3t + 83) 
U, 
It is impossible to identify the dynamics parameters 


215 Dro Pg DO, WD QO; O, O O directly from the record of 
the acceleration of the first story. 


0 The solution is to transform the record from the time 
domain to the frequency domain with a Fourier 
transformation: 


ii,(t) = ¥: [a,cos(a,t) + b, sin(@pt)] = ¥.c, sin(@pt + a) 
t n=l n=) 
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A numerical Fast Fourier transformation (FFT) gives the following two graphs, 
which are the representation of the signal in the frequency domain. 


Ch 


40 


u(t) = dic, sin(a,t + @,,) 30 
n=1 
20 





The required dynamics parameters can be easily determined form these two graphs: 


U(t) = 30sin(27-1.58t+ 30°) + 40 sin(27- 4.43t — 30°) + 20 sin(27- 6.41t + 10°) 


ako Sala2 Qj, = 30 A, = 30° 
f, = 443 Hz A2 = 40 0, = -30° 
f; = 641 Hz Az = 20 0; = 10° 
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The eigenmodes can be entirely determined by doing the same operation for the 2" and 3% story. 


U, (t) M1 M2 3 
Uy(t) > = 51 p SIN(@t+O) + 5 Go p SIN(@yt+A)) + 5 h3 7 SIN(@3t+ 43) 
Us(t) Pr) P39 P33 


i,| 





The sign of each ¢g; is determined by the phase which is also given by a FFT. 
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Free vibration test (with damping) 


The same method is applied. The signal becomes periodic by 


; artificially adding the same model (See Lesson 3) 
Uy 


a 


FFT 


f h h, 


€,, &, é are determined by band width method or by cutting the signal in the frequency domain and 
performing an inverse FFT to obtain free vibration of a SDF system corresponding to the eigenmode. 
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Pulse force excitation 


A hammer is used to apply a short pulse force. 


p(t) u(t) 


t t 


p(t) and the response are non periodic. They 
become periodic by artificially adding the same 
model. The artificial response is then the steady 
state response of the structure loaded by the 
artificial periodic load. The artificial load can be 
considered as the sum of harmonic terms. Each 





of these harmonic terms gives an harmonic Din A 
response. O,=@, — resonance ,D,,;= 5, 
2 D 
Pan(t) = 2. Pp SIN(@,t + a) a fs 
n=l 
Ui, (t) {Din The dynamic parameters can then be 
U, (t) =» 4D>,;sin(a,t+ @,) obtained by performing a FFT. 





68 =| 
U3 (t) art ‘. D3, 


Pn 
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Pan (t) = 2. Pn SIN(@pt + a) 


n=1 


The load obtained by the hammer 
must contain all the frequencies that 
have to be studied. 


in(t)} =F {Do_} sin(apt+,) 


e =I 
U3 (t) art i D3, 


The FFT of the response is similar to 
the graph obtained with an harmonic 
force test (see page 9.5 ) 
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SIGNAL ANALYSIS - FFT 


A signal can be represented in two different ways: in the time domain or in the frequency 
domain. Each domain has its own interest. The transformation between the two domains is 
performed numerically by a FFT or an inverse FFT. 


FFT . | 
et a,(f,) c,(f,) y(t) = > [a,cos(2zf,t) + b, sin(2xf,t)] 
rn or a 
. Palin) An\In) = > c,sin(2nf,t+a,) 
inv. FFT pai 


The signal in the frequency domain is represented by two graphs (a, b,) or (C, @,). 


y Ch Gn N/2 poi 
points 
N points FFT N/2 points A 
Ny 
foo 
mer TT \ 
: Zs fh, <-> In 

0 ie 0 Af Imax 0 Af Imax 


N points FFT 
—<——$_—_—__ > 
+ 
' inv. FFT 
0 7 


y(t) is defined numerically by N points 
during a time T. 


The FFT requires N=2” m= 1,2,3... 


N 
raie = — 
oe 
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N/2 points 


N/2 points WA 


c, and @, are defined by N/2 points. 


1 N _ rate 
uF Imax a oy >. 5 


Ay = 


In practice, (A; , fra.) are chosen 
and (7,rate) are then determined 
such that N= 2”. 


STRUCTURAL DYNAMICS — 5C1840 


¢ displacements ? 
¢ strains ? 


¢ stresses ? 


restriction : linear structures mi+cu+ku= F(t) 









2D continuum beams — Crem > 
¢ Rayleigh’s Fae 
* mass matrix 


¢ pulse and periodic loads 


i 


Corse) — Cor es) — Cea 
¢ eigenvalue problem 
Experimental 
Methods 


¢ mode superposition 


¢ free vibration 
¢ harmonic loads 
¢ time step methods 








free vibration 
damping 
natural frequency 
resonance 
transient response 
steady state response 
explicit time integration 
implicit time integration 


















eigenvalue problem 
eigenmodes 

mode superposition 
mode damping 
Rayleigh damping 
consistent mass matrix 





lumped mass matrix 
time domain 
frequency domain 








